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Abstract 

In this work we perform a detailed numerical analysis of (1+1) dimensional lattice 4 theory. We explore 
the phase diagram of the theory with two different parameterizations. We find that symmetry breaking 
occurs only with a negative mass-squared term in the Hamiltonian. The renormalized mass wir and the field 
renormalization constant Z are calculated from both coordinate space and momentum space propagators 
in the broken symmetry phase. The critical coupling for the phase transition and the critical exponents 
associated with m#, Z and the order parameter are extracted using a finite size scaling analysis of the data 
for several volumes. The scaling behavior of Z has the interesting consequence that (<j> R ) does not scale in 
1 + 1 dimensions. We also calculate the renormalized coupling constant Xr in the broken symmetry phase. 
The ratio X^/wig does not scale and appears to reach a value independent of the bare parameters in the 
critical region in the infinite volume limit. 
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I. INTRODUCTION 



Over the years the 1 + 1 dimensional 4 theories have been used and investigated for many pur- 
poses, including theoretical and algorithmic developments in novel nonperturbative approaches. 
There is a large body of work that deals with the theory in the continuum starting from the mid- 
seventies till now They involve techniques such as Hartree approximation, Gaussian effective 
potential, post Gaussian approximations, random phase approximation and discrete and contin- 
uum light front Hamiltonian. These studies have been done with a positive bare mass-squared 
(m 2 > 0) and diverging contribution to the mass at the lowest non-trivial order of the coupling 
(0(A)) arising from normal ordering was cancelled by a counter term, effectively dropping the di- 
vergent piece. A phase transition to broken symmetry phase was found at strong quartic coupling. 
Critical value for A/w 2 has been found to lie somewhere between 30 and 60. 

Lattice regularization is naturally suited to determine the phase diagram of a quantum field 
theory. There has been an attempt on the lattice ||2j] to extract the critical value of A/m 2 . This 
calculation was performed with a negative bare mass- squared term in the lattice action resulting 
in the broken symmetry phase at small coupling. The negative mass-squared was converted to a 
positive mass-squared in the infinite volume limit by a renormalization performed after the lattice 
data had been extracted and a critical value for A/m 2 was quoted. 

In this paper we investigate the 1+1 dimensional 4 theory on the lattice. To the best of our 
knowledge, there does not exist any detailed study of the critical region of this theory using the 
nonperturbative numerical program of quantum field theories on the lattice. Our aim is to explic- 
itly determine the scaling behavior of the renormalized mass, the renormalized coupling and the 
field renormalization constant including their amplitudes. We also want to investigate the topo- 
logical sector of this theory in the broken phase. In a companion work yfl we have calculated the 
topological charge using the same nonperturbative techniques and have shown its relation to the 
renormalized parameters in the quantum theory. 

In this theory the quartic coupling has the dimension of mass-squared and a 'physical' (relevant 
in the continuum) quantity to calculate is the dimensionless ratio of the renormalized parameters 
Afl/m^. We determine the phase diagram in the two dimensional bare parameter space which 
agrees with the phase diagram of \2]. We have not found a phase transition from the symmetric 
phase to the broken symmetry phase with a positive mass-squared term in the action. The sym- 
metry breaking occurs in our lattice theory only with a negative mass-squared term. We perform a 
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detailed study of the scaling region of the broken symmetry phase and determine the ratio A^/m| 
which appears to be constant in the scaling region irrespective of the bare lattice parameters. The 
vacuum expectation value of the renormalized field (fa) also seems to be constant irrespective of 
the parameters in the scaling region of the 1 + 1 dimensional theory. We have determined these 
ratios using numerical lattice techniques on a variety of lattice sizes. We have estimates for their 
infinite volume values. 

For a nonperturbative approach like the lattice, the notion of perturbative renormalizability is 
to be replaced by existence of critical manifolds and universality classes. The 4 theories are gen- 
erally believed to be in the same universality class as the Ising model. This notion originally came 
from the Renormalization Group and in two Euclidean dimensions is consistent with a conjecture 
based on conformal field theory [4]. In our investigation we determine the critical exponents of 
(0) and iriR independently in (j) 4 theory and find them to be the same as the Ising values. Another 
important ingredient relevant in our analysis is the field renormalization constant Z which appears 
in the two point correlation function. The critical exponent of Z emerging from our FSS analysis 
in 4 theory is found to be consistent with the exponent of susceptibility in the Ising model. 

Although it is not the ultimate goal of our work, our results provide an independent confirma- 
tion of the universality of Ising model and 4 theory (with negative mass-squared) in 1+1 dimen- 
sions using numerical techniques of lattice field theory. However, in an actual calculation, always 
done in a bare theory with an ultraviolet cut-off like the lattice (which also has a finite size), there 
are many important issues still to be resolved, for example, the onset of the scaling region, effects 
of the finite size, possible scaling violation etc. The above needs to be done in each theory for a 
complete understanding of the process of the continuum limit. 

In order to determine the ratio A^/m|, we need to know the field renormalization constant Z and 
the renormalized mass rtiR which can be defined and determined in two ways: 1) the exponential 
fall-off in Euclidean time of the zero- spatial-momentum bare lattice propagators in the coordinate 
space, and 2) the behavior of the momentum space propagators for small four-momenta. We expect 
respective critical exponents corresponding to the renormalized mass and the field renormalization 
constant to agree for the two methods and we verify this in the current paper. In this connection, 
we wish to point out that we have made use of finite size scaling for accurate determination of the 
critical point and verification and determination of the critical exponents. For recent calculations 
in 3+1 dimensional Ising model, see 

Cluster algorithms, known for beating critical slowing down in Ising models, are not directly 
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applicable to the 4 theories. However, owing to a development by Wolff [ 6] using embedded Ising 
variables [7] we have been able to use cluster algorithms in conjunction with the usual Metropolis 
Monte Carlo. 

The plan of the paper is as follows. In Sec. II we define the two parameterizations of the 4 
theory on the lattice followed by section III where we discuss the use of embedded Ising variables 
in 4 theory for use of the cluster algorithms. We show the phase structure of the lattice theory in 
section IV. In section V we present the calculation of the connected scalar propagator in coordinate 
space and momentum space and then in section VI we extract the critical exponents for mass wir 
and field renormalization constant Z and determine the critical coupling using finite size scaling 
analysis. In Sec. VII the renormalized coupling constant Xr and the quantities Xr/iti r and <j)R are 
discussed. Finally in Sec. VIII we conclude with a summary of our results. 

H. <p 4 THEORY ON LATTICE 

In this section we present the lattice action of (1+1) dimensional lattice 4 theory in two differ- 
ent parameterizations. 

A. Parameterization as in the Continuum 

We start with the Lagrangian density in Minkowski space (in usual notations) 



Note that in one space and one time dimensions, the scalar field <j) is dimensionless and the quartic 



(2.1) 



which leads to the Lagrangian density in Euclidean space 




(2.2) 



coupling X has dimension of mass 2 . 



The Euclidean action is 




(2.3) 



Next we put the system on a lattice of spacing a with 




X 



(2.4) 
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Because of the periodicity of the lattice sites in a toroidal lattice, the surface terms will cancel 
among themselves (irrespective of the boundary conditions on fields) enabling us to write 

(d^) 2 = -0^0 (2.5) 

and on the lattice 

^0 = ^2 [0x+ M + 0x- M -20x] • (2.6) 

0x±/i is the field at the neighboring sites in the ±/i direction. Introducing dimensionless lattice 
parameters m 2 , and Ao by m 2 , = m 2 a 2 and Ao = A a 2 we arrive at the lattice action in two Euclidean 
dimensions 

s= -EE****** + P+f ) E E ^ • ( 2 - 7 ) 

X \L L X X 

We shall henceforth call this lattice action the continuum parameterization. 

All dimensionful quantities in the following are expressed naturally in the lattice units, basi- 
cally meaning that they become dimensionless in the lattice formulation by getting multiplied by 
appropriate powers of the lattice spacing a. 



B. Another Parameterization 

A different parameterization in terms of field 4> and parameters K and A, henceforth called the 
lattice parameterization is obtained by setting 

(I> = V2k<Z>, ml = Id, Ao = 6^ (2.8) 

K K L 

where, d = 2 in our case. This leads to the lattice action 

S' = -2k ££4> x 4> x+m + E«5 + *> E(^- ! ) 2 ( 2 - 9 ) 

X /i X X 

where we have ignored an irrelevant constant. 

In the limit A — > °°, configurations with 4> 2 7^ 1 are suppressed. As a result, field variables 
assume only two values 4> x — > ±1 and 5' is reduced to the Ising action Si s i ng with 

Rising = -2k ££4> x 4> x+m . (2.10) 

x /I 
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The lattice action given in Eq. ( 12.91) is invariant under the staggered transformation 

K^-K, 4> x ^4> st , x (2.11) 

where <J> stx = ( — l) x i+ x 2 cj> x . As a result, if K c is a critical point, there exists another criti- 
cal point at — Kc. We have three phases: broken phase for K > K c ((<!>) 7^ 0, (4> s t) = 0), sym- 
metric phase for — K c < K < K c ( (<3>) = (4> st ) =0) and a staggered broken phase for K < 
— Kc ( (4>) = 0, (4>st) 7^ 0). Note that the staggered broken phase is inaccessible in the contin- 
uum parameterization. 

Wherever possible, we have made use of the lattice parameterization to check the implementa- 
tion of the our algorithm since it allows a cross checking of the critical points. 

DDL ALGORITHM FOR UPDATING CONFIGURATIONS 

It is well known that most algorithms become extremely inefficient near criticality (i.e near the 
continuum limit). This phenomenon is known as critical slowing down (CSD). To beat CSD in 
our 4 theory (which has an embedded Ising variable as explained below) we have used a cluster 
algorithm (known to beat CSD in Ising-like systems) to update the embedded Ising variables and 
combined it with the usual Metropolis Monte Carlo algorithm. The variant of the cluster algorithm 
that we have used is due to Wolff [8] and is known as 'single cluster algorithm'. 

To see what an embedded Ising variable is and how it is made use of, note that the part of the 
4 action that responds to a change of sign is 

x,fl 

= -£l0*||k+-Ml <T * <T *+M 

= -^/x^fiffi^/i (3.1) 

x,fl 

where o x = sign(0 x ) is called the embedded Ising variable and J XjX+ ^ = \<j> x \\<l>x+n\ resembles a 
coupling that depends on both position and direction. 

Notwithstanding the resemblance, the above action does not describe an inhomogeneous, 
anisotropic Ising model since the couplings J will vary over configurations. In general one cannot 
update different aspects of a degree of freedom ( for example the modulus and sign of 0) separately 

and independently. Nevertheless, updating the Ising sector (sign of <b) in 4 theory is legitimate 

n 

owing to a result due to Wolff [6] that we will call Wolff's theorem. 
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We describe the theory underlying this procedure, not always easy to find elsewhere. We start 
by explaining Wolff's theorem . 

Consider a group G of transformations T acting on the configurations C of some system: 

C^TC. (3.2) 

Now let us consider the group G as an auxiliary statistical system whose (micro)states are the 
group elements {T}. We define the induced Hamiltonian governing the distribution of the auxiliary 
system by H(TC) where H is just the Hamiltonian of the original system now considered as a 
function of T. 

Let us now define an algorithm W with transition probabilities p(C; T — > T') such that 

1) L T e- H{TC) p(C;T - T') = e - H ( T ' c \ 

2) p{TC J x -> T 2 ) = p{C-T x T -> T 2 T). 

Wolff's theorem states that a legitimate algorithm for updating the original system is 

• Fix C = C\ and T = I (identity transformation), 

• Update 7i — > T 2 using the W algorithm, 

• Assign C 2 = T 2 C\ as the new configuration. 

We shall now demonstrate that if Wolff's theorem is applied to 4 theory with an appropriate 
choice of G, it is equivalent to updating the Ising sector independently. 

Take G = Z% where N is the number of sites. Elements T of can be represented by Ising 
variables. T = {o x } with o x = +1 or — 1 so that 

TC=T{(j) x } = {o x (l) x }. (3.3) 

The induced Hamiltonian is 

H(TC) = -Y,[°x<l>x][ex+n<l>x+n] 

= -Y,\<lhc\\<l>x+p\OxSx<*x+nSx+ll 

x,n 

(3.4) 

where X = \ (j) x \s x and s' x = o x s x . 
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The proof of our proposition follows if we note that the above Hamiltonian is indeed an in- 
homogeneous, anisotropic Ising Hamiltonian and there is a 1-1 mapping between the variables s' x 
and o x . 

We have used Wolff's single cluster variant of the cluster algorithm to update the Ising vari- 
ables. Since the configuration space for the <j)'s is much larger than that for the Ising variables, to 
ensure ergodicity, the algorithm was blended with the standard Metropolis algorithm. The blend- 
ing ratio used was 1:1 i.e every cluster sweep was followed up by a Metropolis sweep. 

We summarize the main steps in the algorithm. We start with some initial configuration for 
the <j) fields. We then update the sign of the fields using Wolff's single cluster algorithm using 
the action (13.11) : We choose some site (seed) at random and select a group of fields (cluster) 
around the seed having the same sign as the field sitting at the seed. The probability for selecting 
a particular field is governed by the action (13 .lb . This process is called growing a cluster. We flip 
the sign of the fields belonging to the cluster (the variables a in (13.11) ) when the cluster is fully 
grown. Finally we execute a Metropolis sweep over the entire lattice updating the full <j) fields. 
This completes one updation cycle. 

Throughout this paper we have used periodic boundary conditions. However, in the companion 
work [3] dealing with topological charge we have used antiperiodic boundary conditions where 
cluster algorithms do not work J^]. 



IV. PHASE STRUCTURE 

The phase structure is determined by looking at the order parameter (0) which takes a nonzero 
value in the spontaneously broken phase. With the cluster algorithm however, since the sign of the 
field of all the members of the cluster are flipped in every updation cycle the algorithm actually 
enforces tunneling between the two degenerate vacua in the broken phase. As a result, as an 
artifact, the average of over configurations, i.e., the expectation value becomes zero. Thus to 
get the correct nonzero value for the condensate we measure where = Vol J ime £ (j) x . To 

sites 

understand the mod let us consider a local order parameter (<j) x ). Since the configurations will be 
selected at random dominantly from the neighborhood of either vacua in the broken phase, (<j) x ) 
will vanish when averaged over configurations thus wiping out the signature of a broken phase. If 
one uses (\(j) x \} as the order parameter then in the broken phase it correctly projects itself onto one 
of the vacua yielding the appropriate non-zero value. The use of this mod, unfortunately, destroys 
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FIG. 1 : Phase diagram for continuum parameterization 



the signal in the symmetric phase completely by wiping out the significant fluctuations in sign. 
However if we choose to use ( | volume £ 0*1)' ^ correctly captures the broken phase as well as the 

sites 

symmetric phase. While the sign fluctuation over configurations are still masked, the fluctuations 
over sites survive producing (|0 1) =0 correctly in the symmetric phase. 

The phase diagram for the continuum parameterization obtained for a 512 2 lattice is presented 
in Fig. [T] Classically, spontaneous symmetry breaking (SSB) occurs for negative m 2 ,. For small 
negative m 2 ,, as two minima are shallow and very close to each other, quantum fluctuations can 
restore the symmetry. So, larger negative values of m 2 , are required for SSB to take place. Conse- 
quently, the phase transition line is found in the negative mk semiplane. Our phase diagram agrees 

nfi n 

qualitatively with that obtained for much smaller lattices in [10, 11]. In [2] the authors have ex- 
trapolated their results to infinite volume. We find that our 512 2 lattice results are as good as the 
infinite volume result in [2] . 

In lattice parameterization, phase diagram obtained for a 100 2 lattice is presented in Fig. Eta). 
We have restricted ourselves only to K > region. The symmetry of the phase diagram for K < 
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FIG. 2: (a) Phase diagram for lattice parameterization, (b) Manifestation of the staggered symmetry of 
lattice parameterization of the action. 

and k > is evident from the behavior of (4>) and (4> st ) as a function of K, shown in Fig. Etb). 

As mentioned in the Introduction, in the continuum version of 1 + 1 dimensional 4 theory 
with a positive mass-squared term in the Hamiltonian, there have been many attempts to calculate 
critical couplings for phase transition from the symmetric phase to the broken phase [1]. We have 
investigated the phase diagram of the lattice theory in the region of positive mass-squared and 
have been unable to detect a phase transition in this region of the parameter space. In Fig. |3]we 
show the measurement of the mass gap mr extracted from coordinate space propagator in positive 
mass-squared region. We find that the mass gap niR monotonically increases with the coupling. 

An alternative way to probe the same region of the parameter space of the continuum parame- 
terization is to perform simulations with the lattice parameterization of the action. From the phase 
diagram for the latter presented in Fig. E£a), we reconfirm the absence of phase transition for 
the lattice theory in the positive mass-squared region in the continuum parameterization since this 
whole region can be mapped onto the symmetric phase in the lattice parameterization using the 
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FIG. 3: Mass gap for + ve wiq. 

transformation Eq. (12.8b . 

V. CALCULATION OF PROPAGATOR 

We have made use of two-point connected correlation function to calculate the fundamental bo- 
son mass and field renormalization constant. We have carried out our simulation both in coordinate 
and momentum space. 

A. Coordinate space 



In coordinate space, 2-point connected correlation function G c (x,xo) is given by 
G c (x,x ) = ((j)(x,t)<j>(xo,to))-((l)(x,t))(<j>(x^,to)) 

= (<Kx,*)<K*Mo)>-(0) 2 - 



(5.1) 
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FIG. 4: niR from coordinate space propagator for dif- FIG. 5: Z from coordinate space propagator for dif- 
ferent lattices (a) 32 2 , (b) 48 2 , (c) 64 2 , (d) 80 2 , (e) ferent lattices (a) 32 2 , (b) 48 2 , (c) 64 2 , (d) 80 2 , (e) 
96 2 and (f) 128 2 (m 2 = -0.5) 96 2 and (f) 128 2 (m 2 = -0.5) 

For the derivation of the last equation translational invariance has been assumed. As explained 
before, we actually calculate instead of (0) in Eq. (15.11) The notation x may be confusing in 
1+1 dimensions; however, it is kept to distinguish between the spatial and temporal directions. 

We have extracted the renormalized scalar mass iur (pole mass) and the field renormalization 
constant Z from 



G r (t) 



e ~m R t _|_ e -m R (L-t) 



+ higher states. (5.2) 



2niR 

where G c (t) is the zero spatial momentum projection of the 2-point connected correlation function. 
The second exponential term in the RHS of the above equation is due to the periodicity of the 
lattice. 

In this calculation, to ensure thermalization, we have discarded the first 10 6 configurations 
before starting our measurements. Measurements were carried out on 100 bins of 2 x 10 5 config- 
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urations. In each bin, to fulfill the requirement that measurements be made on statistically inde- 
pendent configurations, measurements were performed every tenth configuration (hop-length). 

Figs. |U and E] show m R an d Z extracted from the connected propagator in coordinate space for 
lattices of six different sizes as function of Xq for fixed m 2 , = —0.5. The figures clearly demonstrate 
scaling of mr and Z as one moves towards the critical point. The critical point is given roughly 
by the dip in each curve. These dips are clearly not at the same place for mr and Z at the smaller 
lattices. A phenomenological finite size scaling analysis has been done in the next section on these 
data and the data obtained from momentum space propagators to calculate the critical exponents 
and the critical coupling in the infinite volume limit. Around the critical region in Figs. |4|and|5l 
all the curves have a thick appearance because of the proximity of the many data points with the 
associated errors shown for each point. Outside the scaling region, a region not of interest to us, 
the data is relatively sparse and we also have suppressed the error bars for them. Unlike in the 
Ising model (Aq — > °°), at finite Aq the (|0|) takes on large values outside the scaling region in the 
broken phase. In calculating the connected propagator, one performs a subtraction between the 
large expectation values of two quantities measured independently. This enhances the error bars 
outside the scaling region in the broken symmetry phase. 

B. Momentum space 

Connected propagator in momentum space is 

G c (p)=^[(^ )-(^)m. (5.3) 

X 

To improve statistics in numerical simulation, averaging over source y is performed. 

G (p) = ^I^- y) K^>-(^)<0v>] 

x,y 

= (^^x$yCosp(x-y)^-m\ 2 8(p). (5.4) 
At small momenta, the momentum space propagator behaves as 

m R' + P 

where, p 2 = 4^ sin 2 (^j-) with ji = 1, 2 is the dimensionless lattice equivalent of the momentum 
M 

square in the continuum. 
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FIG. 6: (a) in the critical region for different L from momentum space propagator, (b) Z' in the critical 
region for different L from momentum space both for m§ = —0.5 



From the intercept of inverse propagator on the ordinate and slope at p 2 — 0,m R i and Z' can be 
determined. However it is only near the critical coupling that the pole of the propagator is actually 
near zero and it is here that m R i approaches the pole mass m R . A similar argument applies to Z' 
and Z. We have thus calculated the momentum space propagators only near the critical point. 

Since the calculation of momentum space propagators were extremely time consuming, we 
had to restrict ourselves to smaller lattices. For thermalization 10 5 configurations were discarded 
before starting the measurements. The number and size of bins as well as the hop-length were the 
same as that for coordinate space propagators. 

As the inverse propagator was found to be nonlinear in the small momenta region, one could use 
only the lowest few momentum modes for the determination of m R i and Z' . For this calculation 
we took the lowest 2 and 3 modes excluding the zero momentum mode because the connected 
propagator value for the zero mode is prone to relatively large statistical error arising from the 
subtraction, in the critical region, between two quantities measured independently [5]. In Figs. 
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|6] (a) and (b), we have presented m R i and Z' extracted from momentum space propagator with 
Mq = —0.5 for four different lattices. Results obtained using 2-point and 3-point fitting are found 
to be quite close to each other. The 3-point fitting is more stable and has been used in our study of 
finite size scaling in the next section. 



VI. CRITICAL EXPONENTS AND CRITICAL COUPLING FROM FINITE SIZE SCALING 
ANALYSIS 

In this section, we perform a phenomenological Finite Size Scaling (FSS) analysis of our data to 
extract critical exponents and critical coupling. Let us first briefly summarize the main aspects I12I 
of this FSS analysis. In a finite size system, there are, in principle, three length scales involved: 
correlation length <§, size of the system L and the microscopic length a (lattice spacing). FSS 
assumes that close to a critical point, the microscopic length a drops out. According to FSS J3], 
for an observable Pl (whose infinite volume limit displays nonanalyticity at the critical point Aq), 
calculated in a finite size of linear dimension L, 

ft(T)/P»(T)=/(L/£»(T)) (6.1) 

where x = (Aq — Ao) /Aq and the function / (commonly known as scaling function) is universal in 
the sense that it does not depend on the type of the lattice, irrelevant operators etc. It does depend 
on the observable P, the geometry, boundary conditions etc. For fixed L as x — > 0, strictly there is 
no phase transition. Consequently Pl(t) is not singular at Aq = Aq. Near the critical point we have, 
^oo(t) = A^x~ v where v is the critical exponent associated with the correlation length. Suppose, 
near the critical point, Poo(x) = ApX~ p . Then from Eq. (16.11) 

P L (x)=A P x-P f(A^Lr v ). (6.2) 

Since Pl{t) should have smooth behavior as x — > 0, a simple ansatz for / may be taken as 
f(L x v ) ~ (L x v )Pl v so that P L (x) does not blow up as x -> 0. Thus f(x) ~ C.W V as x -> 0. 
Alternatively, we may write 

P L {x) =A P A/ /V LPl v g{A- l/v tL 1 ^), (6.3) 

where g is another scaling function. Since Pl(t) should have no singularity as x — > 0, L finite, we 
have, g(x) — > constant as x — >• 0. 
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FIG. 7: Determination of the critical exponents and critical coupling from finite size scaling analysis of 
data for m.R (or m^) in the critical region. L m R and L m R i are plotted for different L versus Ao with 
(a) rriR extracted from coordinate space propagator data and (b) m R i extracted from a 3-point fit near zero 
momentum of momentum space propagator data 



Thus we have 



LPl v /P L {x) =A P l A P E /V F(A- l/v tL l / v ) 



where the function F is the inverse of the scaling function g and F(A^ ' v t L 1//v ) — > a constant 
as T — > for finite L. So, we can write 



lP' v /P 1 



A- 1 A p/v 



C P + D P A^ 1/V tL'/ v + 

as T — > (Cp and Dp are universal constants in the same sense as the finite size scaling functions). 
The utility of Eq. (16.41) is that if we plot L p / V /Pl (t) versus the coupling Xq for different values of 
L, all the curves will pass through the same point when T = or equivalently Aq = Aq H14I1 . These 
ideas provide us with a very good method for evaluating the critical point and checking the critical 
exponents. 



(6.4) 
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FIG. 8: Determination of critical exponent (v) for mass gap obtained from coordinate space propagator for 
ml = -0.5. 

We have performed finite size scaling analysis for the observables , m R ~ l (or m R r l ) 

and Z _1 (or Z /_1 ). The critical behavior of ((f)), m R and Z may be written as 

<0> = V f 
m R = 1 T v 
Z = A z l t 71 . (6.5) 

From the general expectation that in 1 + 1 dimensions, 4 theory and Ising model belong to the 
same universality class, we have used the Ising values for the corresponding exponents as inputs 
in our FSS analysis. Thus, j8 = 0.125, v = 1 and r] = 0.25. 

According to the discussion following Eq. (16.4ft . we have plotted LPl v /Pl{t) as a function of 
Ao near Aq for different L, with Pl(t) = m^ _1 , Z _1 and respectively. Figs. Ha) and (b) 

show the plots of Lyhr and Lm R < for different values of L against Ao with m R and m R i obtained from 
coordinate and momentum space propagators respectively. All results are with = —0.5. We 
can clearly identify the critical point with remarkable precision from these plots and this agrees 

17 




FIG. 9: Determination of the critical exponents and critical coupling from finite size scaling analysis of 
data for Z (or Z') in the critical region. Z L and Z' L 25 are plotted for different L versus Ao with (a) Z 
extracted from coordinate space propagator and (b) Z' extracted from a 3-point fit near zero momentum of 
the momentum space propagator data. 

well with the value shown in Fig. [T] Since in this case p = v, we cannot verify the critical 



exponent for m R (or m R i) solely from these plots. However, on differentiating [ 14], Eq. (16.41) (with 
Pl(t) = l/m,R) with respect to Ao we get 

-4- (L m R ) = 5L 1/V as x ->• 0, 

where B is a constant. The exponent v can be computed easily by taking the logarithm of the 
above equation: 

log — (L m R ) = logfl + - logL. (6.6) 

OAq V 

Determination of v using Eq. (16.61) is presented in Fig. [8j In this calculation, we have used the lat- 
tice data obtained from coordinate space propagator for four largest lattices (64 2 , 80 2 , 96 2 , 128 2 ). 
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FIG. 10: (|0|) L 0125 for different L versus Xq. 



FIG. 11: Extraction of critical exponent associated 
with the order parameter. 



The value of v extracted from our fitting is 1.01 ±0.18 which is consistent with the Ising value 
(v = 1) to our numerical accuracy. 

Fig- El( a ) and (b) also give the value of the universal constant Q appearing in Eq. (16.4b . 
Applying Eq. (16.41) for the case of itir, we find at x = 



L m,R = A^ 1 A| C| = Q 



(6.7) 



because p = v. Discarding the 32 2 data which seem not to conform to FSS, from Fig. |7](a) we 
find that Q ^5.1. Data in Fig. |7](b) is noisy; however, it still gives a value around 5.5. 

Plots of Z L 25 and Z' L 25 against Ao for lattices of different lengths L with Z and Z' computed 
from coordinate and momentum space propagators are presented in Figs. |^a) and (b) respectively. 
We also present the plots of (\(j>\) L 0125 against Ao for different lattice sizes in Fig. [TO! Critical 
points obtained from all these FSS plots are very close to each other. These plots also provide a 
very good confirmation of the critical exponents for Z and order parameter 
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The critical exponent for ( 1 1 ) is also determined by fitting our data for the largest lattice (512 2 ) 
at m 2 , = —0.5 to the corresponding scaling formula for ( 1 1 ) . As is shown in Fig. \TT\ the fit is very 
satisfactory and the results for the critical exponent and the amplitude are: 

j8 =0.1233 ±0.0007 and A^ 1 = 0.9811 ±0.0032 (512 2 lattice) (6.8) 

Also, the critical exponent obtained from (|0|) data with m 2 , = —1.0, although not shown here, 
agrees well with that extracted for m 2 , = —0.5. 

Within our numerical accuracy, critical exponent for calculated in the two different ways, 
namely, the FSS and the direct fit of the scaling formula on the 512 2 lattice, are very close to each 
other, indicating that 512 2 lattice is as good as the infinite system. 



VII. RATIOS X R /m 2 R AND (<j) R ) 

We choose the following definition J15I1 of Xr, appropriate for broken phase, in terms of the 
renormalized scalar mass rtiR (or m R i) and the renormalized vacuum expectation value (0r), 

^ = 3 7^2 =3Z 7^|- (7J) 
(<M 2 (0) 2 

It does not require any knowledge of four point Green function which is computationally de- 
manding. The renormalized coupling Xr calculated using niR and Z from coordinate space prop- 
agator for six different lattices using the above method is presented in Fig. [T2J Error bars are not 
shown outside the scaling region for the reason explained earlier. As evident from the figure, the 
renormalized coupling is close to the tree level result in the weak coupling limit. However, Xr 
deviates noticeably from the tree level expectation at stronger couplings; it has a scaling behavior 
in the critical region and actually vanishes at the critical point modulo finite size effects. 

In the previous section, we have already shown that the results of our numerical analysis are 
consistent with the Ising values of the critical exponents, namely, /3 = 0.125, V = 1 and 77 = 
0.25. This has the interesting consequence that in 1+1 dimensions, the ratio A^/m| (or Xr/m^,) is 
independent of the bare couplings in the critical region, as follows: 

Xr/mr 2 = 3/(<M 2 = 3Z/(0> 2 (7.2) 

~ m-wim-**) 2 * = w-Ao)°, (7.3) 

using 77 = 2/3 = 0.25. 

20 




FIG. 12: Xr for different L from coordinate space propagator for = —0.5. 

In Figs. [T3]and[l4|for m,Q = —0.5 and = —1.0 respectively, we plot the quantity (fa) with 
Z evaluated from coordinate space propagator data, versus Ao for six different lattice volumes for 
a set of bare couplings close to and including the critical region. These figures are consistent with 
Eq. (17.31) which shows that in the infinite volume limit (fa) is independent of the bare couplings. 
The figures show that for larger lattices the value of (fa) gets close to unity along a plateau region 
just away from the critical point in the broken phase and quickly goes to a value close to zero 
on the symmetric phase side. Judging from the trend in these figures, we expect the curve in the 
infinite volume limit to take the shape of a step function at the critical point with (fa) dropping 
from around unity to zero as it passes the critical point from the broken symmetry phase to the 
symmetric phase. 

One can also try to take the infinite volume limit of (fa) in the scaling region (as will be shown 
in the following for the ratio Xr/tt^ in Figs. fT5landfT6T). From our extrapolations, although not 
shown here and already quite apparent from Figs. [T3landlT4l this value seems to be very close to 
unity. 
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FIG. 13: versus Xq for different L at FIG. 14: versus Ao for different L at 

m2 = -0.5. m 2 = -\.0. 

Eq. (17.21) then immediately tells us that the infinite volume limit of Xr/yh^ in the scaling region 
would be close to 3. This is what is indicated in Figs. [l5]and[T6l The significant error bars in 
our data result mostly from inaccuracies in the determination of the field renormalization constant 
Z and do not permit us to take a more accurate infinite volume limit. However, the trend is quite 
unmistakable. 

We notice that for both of the two Figs. [I3]and[l4j the curves for different volumes meet at the 
same value of around 0.65 of (<j>a) at the critical values of Xq (~ 1.93 for = —0.5 and ~ 4.46 
for iuq = —1.0). To explain this, we need to look at Eq. (16.41) from which one can write, at T = 
and finite volume, 



Factors of the lattice linear dimension L cancel between the numerator and the denominator in 
the above equation. From the scaling laws Eq. (16.51) we find that the ratio y/Az/Aq, is the infinite 
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FIG. 15: X^jm\ versus 1 /L for different Ao at FIG. 16: A#/ra^ versus 1 /L for different Ao at 

ml = -0.5. ml = -1.0. 

volume limit of which we find to be constant irrespective of the parameters of the theory. In 
addition, because the coefficients Ca, and Cz are constants (they may, however, depend on things 
like the boundary conditions etc.), the value of (fa) at T = is also constant irrespective of the 
parameters and volume, as shown in Figs. [13] and [14] for two sets of bare couplings. 



VIII. CONCLUDING DISCUSSION 

Our investigation of the 1+1 dimensional 4 theory on the lattice has certainly turned out to be 
more challenging and absorbing than what one would generally expect for a lower dimensional 
theory. 

On the algorithmic front, since the Metropolis algorithm was extremely inefficient near the crit- 
ical region and for our study we required to obtain a large number of uncorrelated configurations, 
we had to incorporate the cluster algorithm. Cluster algorithms are generally applicable only to 
Ising-type systems. In our case, we used a result due to Wolff [6] to apply the cluster algorithm to 
the embedded Ising variables in the 4 theory. To update the radial modes of the fields the standard 
Metropolis algorithm had to be blended with the cluster algorithm. 

We explored the phase diagram of the lattice theory in two different parameterizations. We have 
found that symmetry breaking occurs only with a negative mass-squared term in the Hamiltonian. 

We needed a large number of configurations to get numerically stable results for the connected 
propagators in the broken phase. Away from the critical point, the magnitude of the field is large, 
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and the connected propagator which would be a relatively small number had to be extracted from 
the subtraction of two large numbers. In addition, the momentum space propagators showed signs 
of curvature for small lattice momenta, a fact which made the determination of the renormalized 
mass and the field renormalization constant a tough one and we had to be as close to the zero 
momentum as possible. 

Using a definition appropriate for broken symmetry phase, we have calculated the renormalized 
coupling Xr which invloves mr, Z and (0). At weak coupling limit our result is close to the tree 
level result but deviates significantly in the strong coupling regime. 

We have used the finite size scaling analysis to determine the critical point and verify and 
ultimately determine the critical exponents associated with m Rl Z and (0). Verification of critical 
exponents for niR and Z are performed using the data for both coordinate and momentum space 
propagators. Apart from verifying the critical exponent for (0) using FSS analysis, we have also 
independently determined this quantity by fitting our data for a large enough lattice (512 2 ). Our 
results are consistent with the expectation that in 1+1 dimensions the 4 theory and the Ising model 
are in the same universality class. 

One of the most important observation in 1+1 dimensions is that the field renormalization con- 
stant scales with a particular critical exponent, something that does not happen in 3+1 dimensions. 
This has the interesting consequence that the renormalized field does not scale and in the infi- 
nite volume limit, it drops from a value approximately around unity to zero abruptly as we pass 
from the broken symmetry phase to the symmetric phase. Moreover, the ratio of the renormal- 
ized quartic coupling to the square of the renormalized mass also does not scale and appears to 
be independent of the bare parameters in the scaling region. Numerically this ratio seems to ap- 
proach a value around 3 in the infinite volume limit. However, our infinite volume extrapolations 
are approximate due to large finite size effects and systematic error in the calculation of the field 
renormalization constant Z. 

For reliable extrapolation of the above amplitude ratios to infinite volume there exist methods 
Il6 ] which we have not tried in the present work. We need to have smaller error bars especially 
on the niR and Z data and this can be taken up in a future work. 
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